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We preform large-scale ab initio density functional theory calculations to study the lattice strain 
and the vibrational properties of colloidal semiconductor core-shell nanoclusters with up to one 
thousand atoms (radii up to 15.6 A). For all the group IV, III-V and II-VI semiconductors studied, 
we hnd that the atom positions of the shell atoms, seem unaffected by the core material. In 
particular, for group IV core-shell clusters the shell material remains unstrained, while the core 
adapts to the large lattice mismatch (compressive or tensile strain). For InAs-InP and CdSe-CdS, 
both the cores and the shells are compressively strained corresponding to pressures up to 20 GPa. 

We show that this compression, which contributes a large blue-shift of the vibrational frequencies, 
is counterbalanced, to some degree, by the undercoordination effect of the near-surface shell, which 
contributes a red-shift to the vibrational modes. These findings lead to a different interpretation 
of the frequency shifts of recent Raman experiments, while they confirm the speculated interface 
nature of the low-frequency shoulder of the high frequency Raman peak. 

PACS numbers: 61.46.Df, 63.22.-m,73.22.-f 


I. INTRODUCTION 

Golloidal semiconductor nanoclusters (NCs) or quan¬ 
tum dots (QDs) have attracted a great deal of attention 
due to their applications in the fields of optoelectron¬ 
ics, spintronics, photovoltaic or bio-labelingIKm and to 
their potential to emerge as the key components of the 
next generation displays^. One important milestone was 
the synthesis of core-shell quantum dots, which contain 
at least t wo sem iconductors arrangement in an onion-like 
geometr> ^^ l ^ l ^^^ ^. In this, the core material is surrounded 
by a different shell material in order to reduce the in¬ 
fluence of the possibly imperfect surface onto the core. 
Indeed, semiconductor core-shell NGs with a hi gh ph oto- 
luminescence quantum efficiency were reported ^ 

The theoretical modeling of these type of struc¬ 
tures, performed at the level of the effective mass 
approximatiorP^, th e k.p method^, the tight- 
binding awroactP^^^ or the empirical pseudopotential 
methocP^^ were done assuming perfect, unstrained, and 
unrelaxed atomic positions. Indeed, colloidal quantum 
dots in their fluid environment were often considered 
as unstrained, in direct contrast to their self-assembled 
(Stransky Krastanov) omologues that are known to only 
exist because of the presence of strain. The estimate of 
the atomic relaxation and ensuing strain is not straight 
forward in colloidal quantum dots as the use of contin¬ 
uum model^^^Jf^ or empirical force fielcP^ would require 
insights about the surface effects. These are difficult 
to quantify as the surface relaxes inwards, shortening 
the bond-length, but these bonds are weakened. Only 
very recently, ab-initio large-scale calculations showed 
that the effect of structural relaxation in colloidal QDs 
is importanfP^K^ and a lack of such relaxation leads. 


for instance, to the appearance of unphysical imaginary 
vibrational frequencies^ and red shifts of vibrational 
mode^. 

In this work, we perform large-scale ab initio density 
functional theory (DFT) calculations to study the struc¬ 
tural and vibrational properties of core-shell Si(core)- 
Ge(shell), the inverted Ge-Si, InAs-InP, and GdSe-GdS 
NGs with radii ranging from 13.5 to 15.6 A. We find 
that (i) The shell dictates the atom positions of the en¬ 
tire core-shell NGs. This is especially true for group IV 
Si-Ge and Ge-Si NGs where we see almost no difference 
between the bond length distribution of the core-shell 
NCs and the pure NC made of only shell material. For 
instance, the Ge core in a core-shell Ge-Si NC is com¬ 
pressed to the bulk Si lattice constant (4% compression 
of the bond length), (ii) Both the core and the shell 
are compressed in the InAs-InP and in the CdSe-CdS 
NCs. So the lattice constants of core and shell materi¬ 
als do not undergo the compromise of one having com¬ 
pressive and the other tensile strain, as one may expect, 
(iii) The bond-length distribution in the NCs goes from 
homogeneous (small scattering) throughout the NCs for 
our group IV NCs to inhomogeneous only at the inter¬ 
face for our group III-Vs, to inhomogeneous in the entire 
shell region in our II-Vis NCs. We speculate that the 
long-range Coulombic interaction in the more ionic II- 
VIs is responsible for the large bond distortions, (iv) 
The frequency shifts we obtain, compared to the bulk 
frequencies, for our NCs can all be traced back to two 
fundamental effects: One being the shift of the modes 
according to strain (given by the Griineisen parameters), 
the second being the red-shift created by the underco¬ 
ordination of the near-surface atoms. Both effects tend 
to work against each other since the NCs are typically 
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TABLE I: Comparison of TO and LO frequencies (in cm“^) 
at the r points for bulk CdS and CdSe between the NLCC 
calculations and the full calculation including the d-electrons 
in the valence. 



rYQNLCC 

—state 

I^qNlcc 

state 

CdS 

247.2 

246.5 

296.5 

299.1 

CdSe 

179.4 

174.8 

204.1 

205.2 


compressively strained and the strain effect leads to a 
blue-shift (with a positive Griineisen parameter). 


II. METHOD 


We first construct an un-relaxed NCs by cutting a 
sphere, centered on an atom for group IVs or a cation 
for III-Vs and II-Vis with point group symmetry, out 
of bulk zinc-blende material. Then, the surface atoms 
with only one nearest neighbor bond are removed and 
the surface dangling bonds are terminated by hydrogen 
atoms for group IV atoms and pseudo-hydrogen atoms 
iL* with a fractional charge of 1/2, 3/4, 5/4, and 3/2 
for group VI, V, III, and II atoms, respectively. These 
atomic positions are relaxed using DFT in the local den¬ 
sity approximation (LDA) and Trouiller-Martin norm- 
conserving pseudopotentials with an energy cutoff of 30 
Ry for CToup IV and III-V clusters and 40 Ry for II-VI 
cluster^. In order to calculate the vibrational eigen- 
modes of Cd-VI clusters with up to one thousand atoms, 
we apply a non-linear core correlation (NLCC) to the Cd 
atoms instead of including d-electrons in the valence. The 
calculated transverse optical (TO) and longitudinal opti¬ 
cal (LO) frequencies of bulk CdS and CdSe at the T point 
using NLCC-LDA and fully including the d-electrons in 
the valence are given in Table. [T[ As shown in this ta¬ 
ble, the phonon frequencies calculated using NLCC are 
in close agreement with the full calculations. 

The structures are relaxed until the forces are less than 
3x10“^ a.u.. The dynamical matrix elements are then 
calculated via finite difference and the vibrational eigen- 
modes and eigenvectors are obtained by solving the dy¬ 
namical matrix,!^ 


^ 1 d^V{R) 

^MjMj dRjdRj 




( 1 ) 


where / and J label the atoms, M are the atomic masses, 
V{R) the potential energy, R the atomic positions, U 
the eigenvectors and uj the vibrational frequencies. We 
use ab initio DFT implemented in the CPMD cod^^ to 
optimize the geometry and to calculate the vibrational 
eigenmodes of the NCs. 

In order to analyze the vibrational eigenmodes in terms 
of the core and the shell (surface) contributions, we cal¬ 


culate the projection coefficients 


a 


V _ 

c(s) — 


^Nc(N,) 


( 2 ) 


where, Vc, Ng and N are the core, shell (surface), and the 
total number of atoms; Xj represents the three compo¬ 
nents belonging to atom / and vibrational mode n from 
the 3V component eigenvectors Uj. 

To compare the vibrational properties of core-shell 
NCs with the phonon density of states (DOS) of their 
corresponding bulk, the phonon DOS of bulk InAs, InP, 
CdSe, and CdS are calculated via ab initio density func¬ 
tional perturbation theor}!^. The computation of bulk 
materials is performed using the ABINIT code packag^^ 
with the same pseudopotentials and the same energy cut¬ 
off as their corresponding NCs, and a Monkhorst-Pack 
/c-point mesh is taken as 8x8x8. 

To compare our results with Raman spectroscopy mea¬ 
surements, we calculate the Raman intensities using a 
phenomenological model proposed by Richter et 
Based on this model, the Raman intensity of nanostruc¬ 
tures is proportional to the projection coefficient of the 
vibrational modes of the nanostructure onto bulk mode s 
with a relaxation of the wave-vector selection rul^^^EH, 


I{uj) oc ^ 


nj 


^ (w-W-)2 + (ro/2)2’ 


(3) 


where, ^ is the projection coefficient of the NC mode 
u on the bulk mode n with wave vector q and Fq is the 
natural Lorentzian linewidth (an empirical parameter in 
our work). The coefficient ^ are summed up to Aq = 
1/{2R). 


III. STRUCTURAL PROPERTIES 

In Fig. we show the relaxed atomic positions of the 
Ini 4 iAsi 4 o-In 228 P 208 H 3 oo core-shell NC with a radius of 
15.6 A. In this figure, the core part (InAs) is shown 
as green and purple spheres while the shell part (InP) 
is plotted as green and tan spheres. The small white 
spheres represent the passivants. The core, shell and in¬ 
terface areas are highlighted. 

To describe the structural properties of the NCs, we 
plot the nearest neighbor distances as a function of their 
distance from the cluster center. The results are shown 
for group IV, III-V and II-VI in Figs.[^[^ and[^ respec¬ 
tively. 


A. Strain: Group IV NCs 

For group IV NCs we compare in Fig. the results for 
the core-shell structure with the results for a pure NC 
of similar size made of shell material. The most striking 
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FIG. 1: (Color online) Relaxed atom positions of the 

In79As68-In242P244H3oo coie-shell nanocrystal. The In, As, 
P, and H* atoms are represented as green, purple, tan, and 
white spheres, respectively. The definition of core and shell is 
following the atomic type (As atoms are core-atoms, P-atoms 
are shell-atoms, in the present case). The interface region is 
defined graphically. 


result is that the bond length distribution of both struc¬ 
tures is nearly identical. In other words, the core-shell 
structure is geometrically very similar to a NC made of 
pure shell material, i.e., the shell dictates the structure 
and the core adapts. 

For the Ge-Si core-shell structure^ Fig.^i,b), the softer 
Ge material (bulk modulus of 77.2 GPa^ is fully com¬ 
pressed, by around 5%, to the lattice constant of Si. The 
bond lengths show a very small deviation from the bulk 
values, even close to the surface, no noticeable deviation 
is observed. This is a special feature of Si, compared to 
the other materials. 

For the Si-Ge core-shell structure^ Fig. [^,d ), the 
rather stiff Si core (bulk modulus of 97.6 GPgP^ is ex¬ 
panded by 2-3 %, in order to almost perfectly match 
the bond length distribution of a pure Ge NC. The Ge 
NC, in contrast to the Si NC, has a significant bond- 
length reduction at the surface (around 5%) and even 
for our largest NC (Fig. [^) the bulk bond-length is not 
recovered at the NC’s center. We refer to the bond- 
lengt h red uction at the surface as the undercoordination 
effecW^^ i.e., the missing atomic partners on the vac¬ 
uum side lead to an inward surface relaxation, as common 
in surface physics/chemistry. 

We further note that the bond-length variation is 



FIG. 2: (Golor online) Bond length distribution as a function 
of their distance to the dot center. For (a) Si 465 H 228 (red 
circles), Gei 47 -Si 3 i 8 H 228 (black squares); (b) Si 705 H 3 oo (red 
circles), Ge 28 i-Si 424 H 3 oo (black squares); (c) Ge 465 H 228 (green 
triangles), Sii 47 -Ge 3 i 8 H 228 (blue diamonds); (d) Ge 705 H 3 oo 
(green triangles), Si 28 i-Ge 424 H 3 oo (blue diamonds). The LDA 
bond lengths of bulk Si and Ge are given as dashed lines and 
dotted dashed lines, respectively. 


rather constant across the interface, marked as light- 
green area in Fig. 


B. Strain: Group III-V (InAs-InP) NCs 

The corresponding comparison for InAs-InP is given 
in Fig. 1^ We now notice a different bond-length dis¬ 
tribution between the InAs-InP core-shell structure and 
the pure InP NCs. The difference is, however, still sur¬ 
prisingly small and only apparent for the smaller NCs 
(Figs.[^,b)). In all cases, the core is heavily conmressed, 
by 5.1%-7.1%. For a radius of 16.5 A (Fig. [^)), the 
bond-length of the core are already “converged” to the 
bond length of the shell. So again, the shell dictates the 
structure and the core adapts. The shell in the core-shell 
InAs-InP NC, or the surface-area in the pure InP NC, 
have a reduced bond-length due to the undercoordination 
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Distance to NC center (A) 


FIG. 3: (Color online) Bond length distribution as a func¬ 
tion of the distance to the dot center for (a) InrgAses- 
Ini46Pi72H228 (black square), In225P24oH228 (red circle); (b) 
In79As68-In242P244H300 (black square), In32iP3i2H3oo (red 
circle); and (c) Ini4iAsi4o-In228P208H3oo (black square), 
In 369 P 348 H 3 oo (red circle). The LDA bond lengths of bulk 
InP and InAs are given as dashed lines and dotted dashed 
lines, respectively. 


FIG. 4: (Color online) Bond length distribution as a func¬ 
tion of the distance to the NC center for (a) CdrgSees- 
Cdi46Si72H228 (black squarc), Cd225S24oH228 (red circle), and 
(b) Cd79Se68-Cd242S244H300 (black square), Cd32iS3i2H3oo 
(red circle). The LDA bond lengths of bulk CdS and CdSe are 
given as dashed lines and dotted dashed lines, respectively. 


effect: the surface layers relax inwards. The bond-length 
variation is significant in the interface region and hence 
much larger than in the group IV materials. 

C. Strain: Group II-VI (CdSe-CdS) NCs 

The results for CdSe-CdS NCs are given in Fig. In 
comparison to the previous results for our group IV and 
our group III-V NCs, we can see a clear progression to¬ 
wards a more scattered data: the core-shell structure 
(black squares) shows a bond-length distribution that sig¬ 
nificantly differs from the corresponding NC made of only 
shell material (CdS, red circles). The core CdSe mate¬ 
rial is still compressed by as much as 2.0% but has larger 
bond-length than the core of the corresponding pure CdS 
NC. In the area of the interface, the bond-length varia¬ 
tion is large, going from nearly CdSe bulk bond-length 


to 1.9% below the CdS bulk bond length (again, the un¬ 
dercoordination effect). 

The comparison of Fig. [^) and b) shows as notice¬ 
able difference, especially in the surface/shell region, de¬ 
spite the fact that the NCs have only a small difference 
in radius. This highlights the fact, that the atomistic 
description leads to a shell-by-shell construction of the 
structure with increasing radius. A small radius increase 
can lead to geometrically and chemically rather different 
structures. The pure CdS NCs have a ratio of Cd to S 
atoms of 225/240 in the smaller structure and 321/312 
in the larger structure. So, Cd poor in the first case and 
Cd rich in the second. This should be kept in mind when 
comparing structures with different radii. 
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IV. VIBRATIONAL PROPERTIES 

The vibrational DOS are shown for the different NCs 
in Figures. ii and Two dominant effects lead to 
the vibrational frequency shifts we observe: (i) Strain 
induced shifts. The magnitude of this shift is quantified 
by the Griineisen parameter. For the materials consid¬ 
ered, these are all positive and between 0.89 and 1.89 (see 
Table 0 - The Griineisen parameters for GdS and GdSe 
TO-modes are results of our DFT calculations, as we 
did not find experimental results in the literature. This 
means, for a compression of the structure, the LO and TO 
frequencies shift to higher frequencies (blue shift), (ii) 
The undercoordination effect. The atoms at the surface 
lacking bonding partners vibrate with a lower frequency. 
This red-shift is strongest for atoms close to the surface 
but does penetrate a few layers inside the NGs. 


TABLE II: Griineisen parameters for the optical modeJ^. 
The resnlts marked with an asterisk * are from onr DFT cal- 
cnlations. 



sP Gpl InlE3 InApI Cd^ 

CdSpl 

yTO 

1.02 0.89 1.24 

1.02 0.89 1.44 

1.06 

1.21 

1.89* 

1.37 (1.34*) 

1 .86* 

1.1 (1.43*) 


A. Vibrations: Group IV NCs 


The vibrational DOS is given for “smaller” (693 atoms) 
and “larger” (1005 atoms) structures in Fig.[^ In case of 
the core-shell NGs (panels (a),(b),(e),(f)), the vibrational 
eigenmodes are projected onto a core and a shell area de¬ 
fined in Fig. ^and shown as black and red lines respec¬ 
tively in Fig. ^ For the pure NGs (panels (c),(d),(g),(h)) 
the eigenmodes are projected onto a core and a surface 
area. 

We make few observations: 

1) The vibrational DOS of the core-shell struc¬ 
tures ((a),(b),(e),(f)) is qualitatively a superposition 
of the vibrational DOS of the pure Si and Ge NGs 
((c),(d),(g),(h)). 

2) Gompared to the bulk optical frequencies (blue 
dashed lines at high frequencies in panels (c),(d),(g),(h)), 
the optical peak is red-shifted. In case of the Si NGs (pan¬ 
els (c),(g)) it is entirely due to the undercoordination ef¬ 
fect, as the structure is basically unstrained (see Fig. 
for the bond-length distribution. For the case of the Ge 
NGs (panels (d),(h)), the undercoordination effect (red- 
shift) is counterbalanced by the compression (blue-shift) 
experienced throughout the crystal (see Fig. [^. The un¬ 
dercoordination effect is stronger and the optical peak is 
slightly red shifted. 

3) The high frequency optical peak (originating from 
Si) of the Ge-Si NGs (panels (a), (e)) is significantly blue 


693-atom structures 


(^) ^® 147 '^' 318*^228 

-Core 
—Shell 


1005-atom structures 
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FIG. 5: (Color online) Vibrational DOS of (a) Ge-Si core¬ 
shell NC with R = 13.5 A, (b) Si-Ge core-shell NC with R 
= 13.8 A, (c) pure Si NC with R = 13.5 A, (d) pure Ge NC 
with R = 13.8 A, (e) Ge-Si core-shell NC with R = 15.2 A, 
(f) Si-Ge core-shell NC with R = 15.5 A, (g) pure Si NC 
with R = 15.2 A, (h) pure Ge NC with R = 15.5 A. All the 
vibrational modes were broadened by 0.8 cm“^. The lower 
frequency blue dashed lines in panels (c), (d), (g), and (h) la¬ 
bel the van Hove singularities for the acoustic branches, which 
corresponds to a maximum in the bulk acoustic phonon DOS. 
The higher frequency blue dashed lines show the optical T- 
point frequencies in bulk (517.0 cm“^ for Si and 300.9 cm“^ 
for Ge). The “surface” of the pure Si(Ge) NCs are defined to 
have the same dimension as the “Shell” in the core-shell NCs. 
Note that the passivant vibrations are far remote at much 
higher frequencies, see Ref. |36|for a detailed description. 


shifted with respect to the optical peak of the Si-Ge NGs 
(panels (b), (f)). This shift originates mainly from the 
fact that the Si core (panels (b)(f)) experiences tensile 
strain of up to 3% compared to its nearly unstrained 
value when it is used as shell material (panels (a) (e)). 
This is mainly a strain effect. 

4) For Si and Ge, the Griineisen parameters of the lon¬ 
gitudinal acoustic (LA), LO, and TO modes are positive 
while those of the transverse acoustic (TA) modes are 
negative. Due to the mixing of transverse and longitudi- 
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nal characters in the confined NCs, the positive and neg¬ 
ative Griineisen parameters of the acoustic modes tend 
to cancel each other out. This leads to the lack of shift 
in the acoustic modes that we highlighted in Fig. 
(panels (c)(d)(g)(h)) by comparing the NC’s results with 
the van Hove singularity of the bulk acoustic mode^ 
marked with blue dashed lines at 155 cm“^ for Si, and 
90 cm“^ for Ge. 


B. Vibrations: Group III-V (InAs-InP) NCs 

In a similar fashion as done for group IV, the vibra¬ 
tional DOS is given for InAsP NGs in Fig. We note 
the following: 

1) In contrast to the group IV NCs, the optical peaks 
are strongly blue-shifted compared to the bulk frequen¬ 
cies. This is true for the pure NCs (panels (d) and (h)) 
and for the core-shell structures (panels (b),(c),(f),(g)). 
This can be traced back to the very strong compression 
given in the core (in the case of the InAs core, more than 
6%) and in the shell (1-2% in the InP shell), as well as 
in the pure NC, as shown in Fig Such a reduction 
of 6% in the bond-length corresponds to extreme pres¬ 
sures of around 20 GPa, for the InAs core. A result that 
highlights the need to consider strain effects in colloidal 
QDs. 

2) The surface modes that appear in the phonon gap 
(between acoustic- and optical-type vibrations) between 
240-280 cm“^ in the InP NCs (panels (d) and (h), see 
Refill] for a detailed description of these new modes) are 
in the same frequency range as the core optical modes in 
the core-shell NCs (panels (b) and (f)). 

3) Some vibrational modes are confined around the in¬ 
terface (panels (c) and (g)), but lie in the same frequency 
range as core and shell modes. 


C. Vibrations: Group II-VI (CdSe-CdS) NCs 

In a similar fashion as done for group IV and III-Vs, 
the vibrational DOS for CdSeS NCs are given in Fig. 
We note the following: 

1) The CdSe core “optical” modes at around 200 cm“^ 
(panels (b), (f)) are blue-shifted compared to the corre¬ 
sponding bulk frequencies. This is due to a compression 
of around 2% of the core (see Fig. [^. 

2) The CdS shell “optical” modes are scattered over a 
large range of frequencies between 250-370 cm“^. This 
is similar to the case of pure CdS NCs (panels (d) and 
(h)) and is directly related to the large variation in bond- 
length approaching the surface of the NCs, as depicted 
in Fig. 

3) Interface modes do exist and are of the persistent- 
typ^, i.e., they are in the frequency ranges of the bulk 
CdSe and CdS modes (with the shifts according to the 
existing strain). In other words, as is the case for the III- 
V NCs discussed previously, no additional modes appear 


(a) 


Bulk InAs 
Bulk InP 




(e) 




693-atom structures 


933-atom structures 



Frequency (crn'O 


FIG. 6: (Color online) (a) and (e) phonon DOS of bulk 
InAs and InP (both (a) and (e) panels show the same re¬ 
sults). Vibrational DOS of (b) InAs-InP core-shell NC with 
R — 13.6 A, (c) projection on interface region (enlarged 3.25 
times), (d) InP NC with R = 13.6 A, (f) InAs-InP core-shell 
NC with R — 14.5 A, (g) projection on interface region (en¬ 
larged 4 times), (h) InP NC with R— 14.5 A. All vibrational 
modes were broadened by 0.8 cm“^. The “Surface” of InP 
NCs is defined to have the same dimension as the “Shell” in 
the corresponding core-shell NCs. 


in the intermediate frequency range between the bulk 
bands. 


V. COMPARISON WITH EXPERIMENT 

To allow a comparison with experiment we have used 
the phenomenological model proposed by Richter et 
given in Eq. (|^, to calculated the Raman intensities 
and show the results for CdSe-CdS core-shell CdygSees- 
Cd 242 Se 244 H 3 QQ NC with R = 14.9 A in Fig. While 
the green lines represent the raw data of the intensities, 
we have also used two different values for the broadening 
Fq: 1.5 cm“^ is shown as red line and a broadening of 
8.0 cm“^ is shown as black line. The experimental results 
taken from Ref. [52] are given in panel (b) of the figure for 
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693-atom structures 


933-atom structures 



Frequency (cm"') 


FIG. 7: (Color online) Same as Fig. but for CdSe-CdS 
core-shell NCs with R — lA.l K [(b)'^d)] and with R — 
14.9A[(f)-(h)]. 


comparison. 

To quantify, as well as possible, the origin of the dif¬ 
ferent vibrations, we have plotted the magnitude of the 
vibrational eigenvector \Xj {R) \ as a function of its radial 
position, i.e., distance to the NC center. So for each atom 
I and eigenmode u we obtain one value for \Xj (i?) |. The 
results for all the atoms on one “shell” around the NC 
center (with same value of R) are averaged and plotted in 
Fig. II The vibrations of peak A and A', as well as F and 
F' are qualitatively very similar and only the vibrations 
A and F are plotted. 

Two comments are due up front: (1) The size of our 
NC with a radius of 14.9 A is significantly smaller than 
the experimental size with R around 30.0 A. (2) Peak 
C (221.4 cm“^) is a surface mode of CdS, as can be seen 
in Fig. and in Ref. [53l which we do not expect to see 
in the experiment because of the much larger size and 
corresponding much smaller surface to volume ratio. 

Without considering peak C, we observe for the larger 
broadening (black line in Fig.|^)) a two peak structure, 
similar to the experimental result. Each of these peaks 
is composed of several peaks, as was also deduced from 


the non-lorentzian line-shapes in the experiment!^. We 
now analyze the peaks subsequently. 

Peaks A (282.5 cm“^) and A' (274.6 cm“^) are vi¬ 
brational modes with optical character of CdS^. With a 
broadening of 8.0 cm“^, these two peaks merge into one 
peak with a frequency of 279 cm“^ (Fig.|^, which corre¬ 
sponds to the peak labeled as CdS LO in the experimen¬ 
tal worlP^. Our combined peak is red-shifted compared 
to the experimental results. From the two effects mainly 
responsible for frequency shifts, strain and undercoordi- 
nation, only the latter effect is relevant since the CdS 
shell is already mainly unstrained in our NC (see Fig||). 
The undercoordination effect becomes weaker with in¬ 
creasing shell size and we expect a blue shift of peaks A 
and A' if we go towards the experimental situation with 
a much thicker shell, until they reach the bulk value of 
around 300 cm“^ (Fig.[^)); in good agreement with the 
experiment. A dependence of the blue shift on the shell 
thickness was also reported recentl}!^. The experimental 
peak “LO 297” in Ref. |52| is therefore a bulk-like, un¬ 
strained, CdS peak (see Fig. 1^ to see the localization in 
the CdS shell) without (largej^confinement effect. 

Peak B (264.3 cm“^) is a vibrational mode localized 
at the interface and in the CdS near-surface shell, with a 
small but non-vanishing surface component (within the 
purple area in Fig. see Ref. |53| for a movie of this 
mode). It corresponds well to the “low-energy shoulder 
(LES)” described in the experimental worlP^ and which 
was observed in NCs with different sizes, shapes, surface 
environments and materi als using R aman and photolu¬ 
minescence measurementThis LE S is oft en 
described as a “surface optical” phonon 
which is in rather good agreement with our identifica¬ 
tion. The fact, that is some surface character leads to 
the expectation that it will have some dependence on 
the type of surface passivation, in good agreement with 
some experiments^. 

Peak D (203.5 cm“^) represents a vibrational mode 
with optical character of the CdSe cor^^ combined with 
CdS near-interface contributions (see Eig. (c)) and 
Ref. [53j). This mode has a blue shift (8.5 cm“^) com¬ 
pared to the bulk LO mode in CdSe. The experimental 
blue shift is somewhat larger (214 cm“^), which can be 
traced back to the fact that a structure with a larger 
shell will experience an even larger compression of the 
core and hence a larger blue shift. 

Peak E (197.7 cm“^) are modes with contributions 
from all regions of the NC: core, interface and shell (see 
Eig. El (d)) and Ref. [53|) . They correspond well to the 
“additional intermediate band or IP mode” seen in the 
experimental work. 

Peak F (190.0 cm“^) and E' (184.0 cm“^). These 
modes are CdS shell modes combined with optical sur¬ 
face breathing-type modes (see Eig. [^(e)) and Ref. l53|) . 
They may correspond to the “SO mode” labeled in the 
experimental paper. 

As a final remark, we note that while the large ma¬ 
jority of theoretical work on core-shell QDs neglect the 
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FIG. 8: (Color online) (a) Calculated Raman spectrum (see 
Eq. for CdSe-CdS core-shell NC with R = 14.9 A . Green 
lines represent the raw data of the intensities, red and black 
lines represent the Raman intensities with a broadening of 
1.5 cm“^ and 8.0 cm“^, respectively, (b) Measured Raman 
spectrum of CdSe-CdS core-shell NC with R around 30.0 A 
taken from Ref. El 


FIG. 9: (Color online) Localization of the vibrational modes 
of the Cd79Se68-Cd242S244H3oo NC with R=14.9 A, where 
we have used the magnitude of the eigenmode as indicator 
(see text). The green highlighted area shows the interface 
region and the violet area the region where passivant atoms 
are present defining a surface region. 


VI. SUMMARY 


effect of str ain altogether, some used simple continuum 
model^i^^ leading to the erroneous conclusion, due to 
the absence of the surface effect, that the shell structure 
experiences tensile strain or missed the existence of cer¬ 
tain Raman active modes with interface or surface char¬ 
acter, such as the LES mode (our mode B). The lack of 
the latter has forced more intricate interpretation of the 
experiment in terms of excitonic effects^, which seem, in 
view of the present ab initio results, at least not neces¬ 
sary to interpret the experimental results. A cheap em¬ 
pirical valence force field description would certainly be 
advantageous, but seem presently out of reach due to the 
complexity of surface and interface effects. A reasonable 
valence force field description of the vibrational property 
of III-V NCs was proposed earlieil^for pure NCs. An ex¬ 
tension towards core-shell structures and group II-VI or 
group IV may be worthwhile, although certainly tedious. 


In summary, we have investigated the structural and 
vibrational properties of colloidal semiconductor core¬ 
shell NCs with up to 1000 atoms via ab initio DFT calcu¬ 
lations. We find that the geometry of Si-Ge and inverted 
Ge-Si clusters is determined by the shell part and there 
is no bond length distortion at the core-shell interface 
of these group IV NCs. Also for our III-V and II-VI 
core-shell NCs, the geometry of the shell is similar to 
that of the corresponding pure NCs made of shell mate¬ 
rial. Accordingly, the vibrational DOS of the shell-type 
vibrations in core-shell NCs is very similar to the vibra¬ 
tional DOS of pure NCs made of shell material. Hence, 
the shell experiences no tensile strain but is rather com¬ 
pressed. This fact could be used to improve continuum 
model descriptions that suggested the opposite situation 
(tensile shell). 

For our III-V and II-VI core-shell NCs we find that 
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the bond-length of the core remains very heavily com¬ 
pressed. We also find an obvious bond length distortion 
at the core-shell interface of our III-V and II-VI core-shell 
NCs. This distortion (i.e., scattering in the bond-length 
distribution) extends beyond the interface region all the 
way to the surface in CdSe-CdS core-shell NCs. We link 
the bond-length distortion to the long range ionic interac¬ 
tion, which is strong in the more ionic II-VI NCJ^. This 
large scattering of the bond-length distribution leads to 
a significant broadening of the vibrational bands, which 
is consequently especially prominent for II-VI NCs. 

We also observe a lack of shift in the acoustic modes 
which we trace back to the mixing of LA and TA modes, 
with their nearly co^ensating positive and negative 
Griineisen parameters^. We obtain a blue- (red-) shift in 
the optical modes of the core parts in Ge-Si (Si-Ge) NGs, 
which can be traced back to the compressed (expanded) 
bond lengths. In addition, we find that the interface 
modes of III-V and II-VI core-shell NCs are just within 
the manifold of the core and the shell modes and overlap 
with the surface modes, which will make them difficult 
to identify experimentally. 

Finally, we compared our results with recent Raman 
experiments and give a new interpretation of the results 


based on the lack of tensile strain and the existence of 
the undercoordination effect. The qualitative picture of 
our frequency shifts can be understood based on two, of¬ 
ten competing, effects of compressive strain (blue-shift 
of optical modes) and undercoordination (red-shift). We 
further show that the often discussed low-energy shoul¬ 
der on the Raman spectra originate from interface vi¬ 
brations with small surface character, in agreement with 
most of the experimental interpretations and in disagree¬ 
ment with earlier theoretical models. 
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